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ABSTRACT 

The core-degenerate (CD) scenario for type la supernovae (SN la) involves the merger 
of the hot core of an asymptotic giant branch (AGB) star and a white dwarf, and 
might contribute a non-negligible fraction of all thermonuclear supernovae. Despite its 
potential interest, very few studies, and based on only crude simplifications, have been 
devoted to investigate this possible scenario, compared with the large efforts invested 
to study some other scenarios. Here we perform the first three-dimensional simulations 
of the merger phase, and find that this process can lead to the formation of a massive 
white dwarf, as required by this scenario. We consider two situations, according to the 
mass of the circumbinary disk formed around the system during the final stages of the 
common envelope phase. If the disk is massive enough, the stars merge on a highly 
eccentric orbit. Otherwise, the merger occurs after the circumbinary disk has been 
ejected and gravitational wave radiation has brought the stars close enough for the 
secondary to overflow its Roche lobe radius. Not surprisingly, the overall characteristics 
of the merger remnants are similar to those found for the double-degenerate (DD) 
scenario, independently of the very different core temperature and of the orbits of the 
merging stars. They consist of a central massive white dwarf, surrounded by a hot, 
rapidly rotating corona and a thick debris region. 

Key words: binaries: close — hydrodynamics — supernovae: general — stars: white 
dwarfs — stars: AGB and post-AGB 


1 INTRODUCTION 

Thermonuclear, or Type la supernovae (SNe la), originate 
from the ex plosion of carbon-oxygen white dwarfs — see, 
for in stance iHillebrandt et all (l2013l ) and iRuiz-Lapuentd 
ll2014h for recent reviews on the explosion mechanisms and 
the different evolutionary channels — and are among the 
most luminous events in the universe. Moreover, since they 
arise from the detonation of white dwarfs with a relatively 
narrow range of masses, SNe la typically have similar in¬ 
trinsic luminosities that can be conveniently standardized 
using the properties of the light curves, and consequently 
play a primary role in cosmology. However, the precise na¬ 
ture of the progenitor systems that give rise to SNe la re¬ 
mains still to be unveiled, and is a topic of active research. 
The fact that sev eral peculiar SNe la have been discov- 
ered so far (e.g., iBildken et al.l l2007l : iPerets et al.l l2010l : 
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[Jordan et al.ll2of^ : iFolev et al.ll2013l ) suggests that differ¬ 
ent progenitors and/or explosion mechanisms could trigger 
such outbursts, resulting in normal and peculiar SNe la. For 
instance, the variety of peculiar SNe la and their range of 
properties, has led to explor e sub-Chandrasekhar mass ex ¬ 
plos ions — see, for instan ce, Wang. Justham fc HarJ (l2013l) 
and lScalzo, Ruiter fc SimI ll20l4) . 

It has long been suggested that a white dwarf in a binary 
system — with either another white dwarf — through the 
so-called double-degen erate (DD) scenario llWebbinkI 19841 : 
llben fc Tutukovlll984l) — or a main-sequence or red giant 
companion — through the single-degenerate (SD) scenario 
ilWhelan fc Ibeniri973l: lNomotall982l: Ih an fc Podsiadlowskil 

l2004h — could give rise to a SN la event. More r ecent 
scenarios — see the reviews of Wang fc HarJ (|2012|) and 
iMaoz, Mannucci fc Nelemanj ll2014l) — are variants of the 
two previously mentioned ones, and include the double det¬ 
onation scenario, where helium-rich material is accreted by 
the carbon-oxygen white dwarf fe.g.. lShen fc Bildstenll2()()9l : 
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[Ruiter_et_aJjl2011 ), t he double white dwarf collision scenario 
_|Rasld£_et^Jj20M; iRosswog et al.l l2009l : iThomps^ I2OIII : 
iKushnir et alJl2013f). and the core-deRenerate (CD) scenario 
dSparks fc Stecherjl 19741: Livio &: RiessI 20 0^: Kashi fc Sokerl 


I2OIII: lllkov fc Sokerl 120121 , |2013|: ISoker et alJ l2013h to be 

studied here. Within the CD scenario the white dwarf 
merges with the hot core of a massive asymptotic gi- 
ant br anch (AGB) sta r (iLivio fc Riessl[2003l : iKashi fc Sokerl 
2OIII: lllkov fc Sokerl I201I |2013|) . In a recent paper 


Tsebrenko fc Soker~(l2015bll argue that approximately 20%, 


and likely man y more, of all SNe la c ome from the CD sce - 


disk is not massive enough, the merger takes place after the 
circumbinary disk has been already ejected. In this case the 
driver of the merger is again the emission of gravitational 
waves, so the orbits will have a small eccentricity and the 
temperature of the core of the post-AGB star will be low. 
In snmmary, in the CD scenario a considerably large range 
of eccentricities of the binary system and of temperatures of 
the naked core of the AGB star are expected, at odds with 
what it is foreseen for the DD channel. 

While several simulations of the DD progenitor channel 
have been alr eady pe rformed over the last two decades 


l [Benz et al.l 1990l: ISegretain. C hab rier fc Mochkovitchl 


core of a n AGB star and a white dwarf can indeed be very 
common (iBriggs et al.l2oi5l ~). A table summarizing these sce¬ 
narios with their main advantages and drawbacks can be 
found in iTsebrenko fc Sokerl (l2015bl l. To these we add a note 
on the recently pro posed singly-evolved star (SES) scenario 
llChiosi et al.ll^M^ l. 

The DD and the CD scenarios share many similari¬ 
ties, and in particular they have a common initial evolu¬ 
tionary path. Typically, in a low- or intermediate-mass bi¬ 
nary system the more massive (primary) component fills its 
Roche lobe while crossing the Hertzsprung gap in the color- 
magnitude diagram and the mass transfer episode is stable. 
The primary star then evolves into a carbon-oxygen white 
dwarf. Later, the initially secondary star fills its Roche lobe 
when it becomes an evolved giant star. This time the com¬ 
bination of a deep convective envelope and an extreme mass 
ratio causes an unstable mass transfer episode from the gi¬ 
ant to the white dwarf and a common envelope is formed. 
In the DD channel the common envelope is entirely ejected 
and a final binary system made of two white dwarfs is ob¬ 
tained, where the s econd white dwa rf formed is normally the 
more massive one llWebbinklll984l i. Gravitational wave ra¬ 
diation finally brings both white dwarfs closer and a merger 
occurs. In the CD channel the merger occurs shortly after 
the common envelope phase. 

Note, however, that there are two significant differences 
between the DD and the CD scenarios. In a DD merger both 
white dwarfs are brought together by gravitational wave ra¬ 
diation, a process which lasts for a long time. Consequently, 
it is expected that both components of the binary system 
are cool, and have nearly circular orbits. Moreover, it is also 
foreseen that both white dwarfs will be synchronized. In con¬ 
trast, in the CD scenario the merger of both stars is triggered 
by the interaction with the circumbinary disk. In particular, 
in the CD scenario during the final stages of the common 
envelope phase, ~ 10% of the common envelope remains 
bound to the binary system formed by the AGB star and the 
companion white dw arf, forming a circumbinary disk — see 
iKashi fc Sokerl (1201 ifl . and references therein. Gonsequently, 
it turns out that in a sizable number of cases as the binary 
system transfers angular momentum to the circumbinary 
disk, the separation of the p air decreases and the ecc entric- 
ity of the system increases (lArtvmo wicz et ^ 199ll') whil e 
the core of the AGB star is still hot ( Kashi fc Sokerl[ 20 Ilf) . 


using Monte Carlo binary population synthesis codes con- 

I 997 I: iGuerrero. Garci'a-Berro & IsernI 

2004; 

eluded that the fraction is only 2-10%, but another recent 

Yoon. Podsiadlowski fc Rosswod 

2007; 

population synthesis study suggests that the merger of the 

Loren-Aguilar, Isern fc Garcia-Berrd 20091: IPakmor 

et al. 


2OI0I : iDan et al.ri'20'lll : [Raskin et alJ l2014l i. no full three- 


favoring a merger at the end of the common envelope phase 
or a short time after, during the planetary nebula phase. 
However, it might be also possible that, if the circumbinary 


dimensional simulation of the GD scenario has been done 
so far. The present paper aims at filling this gap, and 
it is organized as follows. In Section [5] we discuss the 
initial conditions considered to simulate the CD scenario. 
It is followed by Section |31 where we briefly explain the 
numerical techniques used in this paper. In Section [T] we 
describe in depth the results of our simulations. Finally, in 
Section [S] we summarize our main findings, and we draw 
our conclusions. 


2 INITIAL CONDITIONS 

2.1 The final phase of the common envelope 

IKashi fc Sokerl (1201 il l estimated that during the last stages 
of the common envelope phase, between ~ 1% and 10% of 
the ejected envelope remains bound to the binary system 
composed of the cold white dwarf and the hot core of the 
AGB star. They furthermore suggested that, due to angular 
momentum conservation and further interaction of the re¬ 
maining common envelope with the binary system, the fall¬ 
back material forms a circumbinary disk around the binary. 
This newly formed disk would later interact with the binary 
system, thus reducing the orbital separation of the pair much 
more than what it was reduced during the dynamical phase 
where the common envelope is ejected, while the eccentric- 
it y of the binary system would increase. Using t he results 
of lArtvmowicz et al.l (ll991^ . IKashi fc Sokerl (1201 il l obtained 
the rate at which the semimajor axis decreases and the rate 
at which the eccentricity of the orbit increases. In addition, 
they also took into account the rate at which the circumbi- 
nar y disk losses ma s s dur ing the interaction. 

IKashi fc Sokerl (I2OIIII found that the final values of the 
orbital parameters critically depend on the amount of mat¬ 
ter ejected from the system through disk winds. Specifically, 
they found that for disks with masses Mdisk 0.12(Mcore + 
Mwo) — where Mdisk, Mcore, and Mwd are, respectively, 
the masses of the disk, the core of the AGB star and the 
white dwarf — the circumbinary disk-binary interaction 
leads to a merger if the speed of the disk wi nd is the escape 
veloci ty. For the particular case studied by IKashi fc Sokerl 
(l201lll . where Mwd = 0.6 Mq, Mcore = 0.77 Mq and 
Menv = 4.23 Mq, the mass of the disk amounts to 4% 
of the mass of the envelope. They also found that for less 
massive circumbinary disks the separation is not reduced 
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Figure 1. Evolution of the orbital parameters of the binary system formed by an AGB star and a white dwarf during its interaction 
with the circumbinary disk, once a large fraction of the envelope of the AGB star has been ejected. The left panel shows the case in 
which a massive circumbinary disk with guigt = 0.12 is adopted, whereas the right panel the case in which Qdisk = 0.10 is considered — 
see text for details. The solid lines and left scale show the evolution of the semimajor axis (a), and the dashed lines and the right scale 
display that of the eccentricity (e). The hollow, large circles indicate the point at which the minimum orbital distance is such that the 
secondary would overflow its Roche lobe for the equivalent circular orbit, = a(l — e) = uecOi whereas the blue solid symbols show 
the orbital distance a when the secondary white dwarf fills its Roche-lobe radius in our simulations, r = tq. Finally, for the right panel, 
the point at which the semi-major axis of the binary system equals the distance at which the secondary reaches the Roche-lobe radius 
for the equivalent circular orbit, a = ueco is indicated by a triangle. Also in this panel the time at which the circumbinary disk has 
been totally ejected is indicated by the solid black circles. 


enough to result in a merger during the interaction. In these 
cases, gravitational wave radiation will reduce the orbital 
separation, a, and the eccentricity, e, further. Consequently, 
the final merger would occur when the orbit of the binary 
system is almost circular. Hence, in these cases it turns out 
that the orbital separation, a, and the distance at closest 
approach, rmin = u(l — e), are very similar. The first mass 
transfer episode occurs when rmin < ro but a > ro, being 
ro the distance between the two stars when the secondary 
overflows its Roche lobe radius in our numerical simulations. 
It is convenient to define here an equivalent circular orbit, 
which is the circular orbit for which R 2 — Rl, where R 2 is 
the radius of the secondary star, a nd Rl is computed using 
the expression of lEggletorJ lll983l) . We call the semi-major 
axis of this orbit oeco- Moreover, since the timescale for 
gravitational wave emission is long, the core of the AGB star 
would be cold dKashi &: Sokerj|201ll). Note, however, t hat al¬ 
though the procedure detailed in Kashi fc Sokerl il201ll ) takes 
into account the tidal interaction during the common enve¬ 
lope phase, to compute the orbital parameters of our system 
at which our SPH simulations start (see below) we did not 
take into account tidal effects once the circumbinary disk 
has been totally ejected from the system. 


2.2 Orbital parameters 

We compute the evolution of the orbital parameters of the 
binary system before the merger occurs following closely the 


procedure of iKashi fc Soke r 1 I 2 OIII ). In parti cular, we con¬ 
sider the same binary as lxashi fc Sokerj ll201ll ) and we adopt 
two different disk masses. Writing 


Aldisk — ^/disk(Afcore “t“ AIwd), 


( 1 ) 


we adopt as reference cases qdisk = 0.12 — a disk with a 
mass just over the critical one — and qaisk = 0.10 — a 
low-mass one. When qaisk = 0.12 the merger occurs just 
before the disk is totally ejected from the system. Hence, the 
merger would occur when the orbit of the binary system is 
still highly eccentric, and moreover the core of the AGB star 
would still be hot, while when gaisk = 0.10 the merger occurs 
long after the circumbinary disk has been totally ejected 
from the system, and hence the orbit will be nearly circular, 
and the core of the AGB ill be cold. 

The time evolution of the orbital parameters for Qdisk = 
0.12 is displayed in the left panel of Fig. [T1 where we show 
the evolution of both the orbital separation (solid line, left 
scale), and of the eccentricity (dashed line, right scale). The 
large, hollow circles correspond to the time instant at which 
the periastron distance is such that the radius of the sec¬ 
ondary would equal the Roche-lobe radius (R 2 = Rl) for the 
equivalent circular orbit, while the solid, blue circles show 
the point at which the mass transfer episode begins — see 
below for details about the mass transfer episode. Note that 
both times are almost coincidential. For this disk, the merger 
occurs at t ~ 270 days after the common envelope is ejected. 

The equivalent diagram for the case in which Qdisk = 
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0.10 is shown in the right panel of Fig. [T] In this case, the 
disk has been totally ejected at early times (the time at 
which this occurs is represented using black, solid circles in 
this figure). For a relatively long time interval the orbital 
separation and the eccentricity of the pair remain almost 
constant, but as time goes on the emission of gravitational 
waves progressively decreases the orbital distance and the 
eccentricity of the pair, being faster for times just before the 
merger begins. Eventually, at approximately fix 10® years af¬ 
ter the common envelope of the system is ejected both stars 
merge. As in the previous case we also show in this panel 
the points where R 2 ~ Rl for the equivalent circular orbit, 
and the points at which the merger starts. Additionally, in 
this panel we also show using red triangles the times for 
which the semi-major axis equals the distance at which the 
secondary reaches the Roche lobe radius for the equivalent 
circular orbit (o = oeco). Note that in this case the eccen¬ 
tricity is small (e ~ 0.1) and that we expect that, given the 
very long timescale of gravitational wave emission, the post- 
AGB star will have a cold core. In cases where merger takes 
place within several times 10® yr, the SN la ejecta might 
catch up with the ejected common envelope, that once was 
a plan etary nebula. This “SN la Inside a PN” is termed a 
SNIP (iTsebrenko fc Sokeill2015dyj . 


2.3 The temperature of the AGB core 

For the sake of completeness for each of the cases described 
earlier we perform two simulations. For the first of these sim¬ 
ulations we adopt a low temperature for the isothermal core 
of the AGB star {T = 10® K), while for the second we adopt 
a hot core (T = 10® K). Admittedly, these are first order 
approximations, since realistic temperature profiles should 
be adopted. However, given that the degeneracy of typi¬ 
cal AGB cores is high we expect that the influence of the 
adopted temperature profiles on the evolution of the merger 
will be small. The white dwarf is always cold, and for it we 
adopt T = 10® K. The adopted temperature has an effect on 
the initial configuration of the core of the AGB star (of mass 
Mcoie = 0.77 Mq). The radii are Rcoie — 1.02 x 10~^ Rq 
and Rcotb — 9.33 x 10“® Rq, and the central densities are 
p™"'® ~ fi.83 X 10® g cm-® and ~ 7.07 x 10® g cm”®, 
for the hot and cold cores, respectively. Note that these dif¬ 
ferences are small, of the order of a few percents, and thus 
given the approximations adopted in this work, have lim¬ 
ited effects on the dynamical evolution of the merger, and 
on the question of whether there is an i gnition upon merger 
dVoon, Podsiadlowski fc Rosswo^l2007l '). 


3 NUMERICAL SETUP 

Both the white dwarf and the core of the AGB star are made 
of equal mass fractions of carbon and oxygen, with fiat pro¬ 
files. We would like to mention here that the outer helium 
layer of ~ 0.01 Mq surrounding the cores of the white dwarf 
and of the AGB star could be potentially important for the 
evolution of the merger. However, with the number of par¬ 
ticles adopted in our SPH simulations this helium shell it is 
practically unresolvable and difficult to follow its evolution 
numerically. 


The SPH code used in this work is an opti- 
mized and updated version of the code used by 
[Guerrero . Garcfa-Berro fc Isernj ll2004ll and later 




Lore n-Aguilar. Isern fc Garcfa-Berrd ll2009lj . 

I^en- Aguilar. Isern fc Garcfa-Berro lioid ') and 

Aznar-Siguan ^ al.l I 2013l j to simulate double-degenerate 


mergers. This code has been extensively described in those 
works. Consequently, here we will only summarize its more 
important characteristics. 

Our SPH code is fully parallel, and allows to run 
simulations with a large number of particles in a rea- 
sonable time. We use the sta ndard cubic spline kernel of 
iMonaghan &: Lattanzid (Il985h . The gra vitational interac¬ 
tion is also softened using this kernel dHernauist &: Kat3 
Il989ll . The search of neighbors and the evaluation of 
the g ravitational fo r ces are performed using an oc- 
tree ( Barnes fc Hud 198^1. The iteration method of 
[Price fc MonaghanI J200'if r is used to compute the smooth¬ 
ing lengths. To deal with shocks we employ a prescrip- 
tion for the arti ficial viscosity based in Riemann solvers 
Jm onaghajj[l997|), and the variable viscosity parameters of 
[Morris fc MonaghanI |l993) are employed. Additionally, to 
suppress artifici al viscos i ty forc es in pure shear flows the vis¬ 
cosity s witch ofjBalsaral (Il99f ) is a lso used. As explained in 
detail in [Aznar-Siguan et al. 1 2013tl . we follow the evolution 
of both the internal energy and of the temperature to en¬ 
sure a reliable evolution of the temperature in a degenerate 
gas. The thermodynamic properties of m atter are computed 
using the Helmholtz equation of state jTimmes fc Swest^ 
[2000lj . Finally, the nuclear network adopted here incorpo¬ 
rates 14 nuclei. The reactions considered in this work are 
captures of a particles, and the associated back reactions, 
the fusion of two carbon nuclei and the reaction between car¬ 
bon and oxygen nuclei. All the thermonuclear reaction rates 
were taken from the REACLIB data base ([Cvburt et al.l 
[ 2 OIOII . Neutrino losses were also include d, being the corre¬ 
sponding rates those of [ltoh et al.l ([l996tl . 

As mentioned before, the first mass transfer episode be¬ 
gins when Tmin < ro < a. To obtain reliable configurations 
when both stars are at closest approach we proceeded as fol¬ 
lows. We first relaxed independent configurations for each of 
the stars of the binary system and we placed them at the 
apastron in a counterclockwise orbit. We then evolved the 
system in the corotating frame to obtain equilibrium config¬ 
urations at this distance. Once this relaxation process was 
finished, we started the simulations in the inertial reference 
frame. In a first set of preliminary simulations we explored 
at which distance the mass transfer episode ensues. This is 
done employing a reduced number of particles (~ 2 x 10^ for 
each star). Once we know which are the orbital separation 
and the eccentricity of the orbit for which we obtain a Roche- 
lobe overflow we computed a second set of simulations with 
enhanced resolution. Since for highly eccentric orbits both 
stars are initially separated by very large distance, to save 
computing time in this second set of simulations we first 
followed the evolution of the system employing the reduced 
resolution and, once the stars had completed one quarter 
of the orbit, these low-resolution simulations were stopped 
and both stars were remapped using a large number of SPH 
particles (a factor of 10 larger). We then resumed the sim¬ 
ulations with this enhanced resolution. This was done be¬ 
cause for systems with large eccentricities the computing 
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Figure 2. Evolution of the binary system for selected stages of the merger, as described in the main text, in the equatorial (left) and 
meridional (right) planes, for the case in which a hot core of the AGB star is adopted and for ^disk = 0.12. The color scale shows the 
logarithm of the temperature, whereas the x, y and 2 axes are in units of 0.1 Rq- The arrows show the direction of the velocity field, but 
not its magnitude, in the corresponding plane. Each panel is labelled with the corresponding time in units of the initial binary period, 
To = 10,126 s. The white dwar f is the obj ect located initially at the right, and later is destroyed. These figures have been done using 
the visualization tool SPLASH jPric 3 I 2007 I 7 See the online edition of the journal for a color version of this figure. 


load required to follow an uninteresting phase of the orbital 
evolution of the pair is exceedingly large. 


4 RESULTS 

As discussed earlier, in this paper we simulate four differ¬ 
ent configurations for the binary system. In all cases the 
mass of the core of the AGB star and the white dwarf are, 
respectively, 0.77 Mq and 0.60 Mq. The first configuration 
corresponds to a binary system that has a disk with a mass 
ratio equal to the critical one of gaisk = 0.12 (see equation 
[TJ, while for the second one we have chosen a disk with a 
mass ratio of guisk = 0.1, which is below the critical value. 
For each one of these cases we have studied the coalescence 
when the temperature of the core of the AGB star is high 
(T = 10® K) and low (T = 10® K). 

For the highly eccentric orbit that results when gaisk = 
0.12 is adopted, the first mass transfer from the white dwarf 


to the core of the AGB star occurs when the distance be¬ 
tween both components of the binary system at closest ap¬ 
proach is Tmin = O.OOaECO, where flECO has been defined 
previously in section [24] as the radius of a circular orbit for 
which the white dwarf overflows its Roche lobe. In comput- 
in g this orb i tal d istance the classical analytical expresion 
of lEggletonl (Il983l i for the Roche lobe radius i?L was used. 
Note that this expression does not take into account tidal 
deformations, and is only valid for circular orbits, while the 
orbits of the binary systems studied here are elliptical and 
our stars are tidally distorted. At this time the semima¬ 
jor axis is a = 1.12 7 ?q and the eccentricity is e = 0.97. 
For the less massive disk with guisk = 0.10 the merger be¬ 
gins when a = 0.038 Rq and the eccentricity of the orbit 
is small e = 0.095. Specifically, the merger occurs when 
rmin = 0.98aECO. Fig [T] illustrates all this. Particularly rel¬ 
evant for the discussion of our results are the values of the 
respective initial periods. For the case of the system with 
gdisk = 0.12 the initial orbital period is To ~ 10,126 s, while 
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Figure 3. Same as Fig. [2] but for guisk = 0.10. In this case the initial orbital period is Tq = 65 s. See the online edition of the journal 
for a color version of this figure. 


for the case in which guisk = 0.10 is used the initial period 
is significantly shorter, To ~ 65 s. 

4.1 Evolution of the merger 

Fig. [5] shows the time evolution of the binary system during 
the merger for the case in which the temperature of the AGB 
core is 10® K and for gaisk = 0.12. Shown are the positions of 
the SPH particles and their color-coded temperatures. The 
left column shows the evolution in the orbital plane and the 
right column in the meridional plane, respectively. In Fig. [3] 
we present the results in the same setting for the case in 
which gdisk = 0.10 is adopted. The white dwarf is the object 
located initially at the right, and later is destroyed, while 
the core of the AGB star is initially located at the left. 

The top panels of Figs. [2] and [3] are at the initial stages 
of the merger, at times t/To — 0.501 for qaisk = 0.12, and 
t/To — 0.640 for gaisk = 0.10, respectively. We remind that 
the orbital periods are To = 10,126 s for gaisk = 0.12, and 
To = 65 s for gaisk = 0.10, respectively. As can be seen 
the white dwarf is significantly deformed, due to tidal in¬ 
teractions, the deformation being larger for the case of an 
eccentric orbit (Fig.[2|. For the second row of panels we have 


chosen slightly larger times, t/To — 0.502 and t/To — 0.74, 
respectively, and show that tidal deformations close to pe- 
riastron are much larger for eccentric mergers. This leads 
to a faster heating of the external layers of the secondary 
star. The third rows show the systems when the system has 
evolved through several passages trough the periastron. For 
gdisk = 0.12 (Fig. [5]) we have chosen the sixth mass transfer 
episode and for guisk = 0.10 (Fig. [3]) the twenty-third one. 
These correspond to times t/To — 1.599 and t/To ~ 22.2, 
respectively. As can be seen, at these evolutionary stages in 
both cases the secondary (the white dwarf) has increased its 
temperature, and more mass has been accumulated on the 
primary core. 

The fourth rows of Figs. [2] and [3l at times t/To — 1.788 
and t/To — 36.5 respectively, display the situation during 
the last orbit, just before the secondary white dwarfs are 
completely disrupted by the cores of the AGB stars. At this 
point, for gaisk = 0.12 the white dwarf still keeps « 80% 
of its initial mass, while for g = 0.10 this percentage is 
somewhat larger, « 90%. The fifth rows correspond to a 
time close to that at which the peak temperature is reached, 
which happens during the final infall of the secondary star 
onto the core of the AGB star. Specifically, we have chosen 
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Figure 4. Trajectories of the centers of mass of the white dwarfs (left panel) and cores of the AGB stars (right panel) for the case in 
which (/disk = 0.12 is adopted (top panels) and that in which (Jdisk = 0.10 is employed (bottom panels). The blue dashed lines correspond 
to the case in which a cold AGB core is adopted, and the red solid lines to that in which a hot core is assumed. The solid black circles 
indicate the initial positions of both stars for each of the panels. 


t/Tf) ~ 1.801 and t/To ~ 38.1, respectively. As it will be 
discussed in more detail later, for guisk = 0.12 these temper¬ 
atures are higher than for Qdisk = 0.10. The reason for this is 
twofold. First, the material of the white dwarf is hotter for 
5 disk = 0.12, a consequence of the previous evolution. The 
second reason is that in this case the matter of the disrupted 
white dwarf falls onto the primary star with a much larger 
radial velocity than for gdisk = 0.10. Consequently, matter 
is compressed more violently in this case. Finally, the last 
rows show the systems at late times. We have chosen times 
t/To — 1.802 and t/To — 38.5, respectively. As can be seen, 


the material of the completely disrupted secondary white 
dwarf is spiralling around the primary star. 

A general feature of our simulations is that the sec¬ 
ondary star, the white dwarf, is totally disrupted after sev¬ 
eral periastron passages. In each of these periastron pas¬ 
sages mass is transferred in a relatively gentle way. Accord¬ 
ingly, in all cases the density and temperature conditions 
for a detonation to develop are not met. As can be seen in 
Figs. [2] and [3] in none of the cases the peak temperature 
reaches 10® K, and hence, although nuclear reactions play 
a role, a powerful explosion does not occur. Nevertheless, 
the helium shells of the white dwarf and the core of the 
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Figure 5. Mass accreted by the AGB core as a function of time, expressed in terms of the initial orbital period. The left panel corresponds 
to (jdiak = 0.12, while the right displays the case in which (/disk = 0.10. Each step in these panels corresponds to a periastron passage. 
The total time span is ~ 25, 315 s for the left panel, and ~ 2, 925 s for the right one. 


AGB star can play a role in inducing a detonation. How¬ 
ever, n one of the studies of merging white dwarfs performed 
so fa r ( Guilloclion_et_^ 2010l: Pibwnslev. Moore fc BildstCT 
2012 ; Moore. Townslev fc BildstenI 1201.31 : IPakmor et al. 


2013 ; IShen fc Moorell2014h has been able to self-consistently 


evolve a binary system with stable helium layers, and con¬ 
sequently this is a topic that, although interesting, de¬ 
serves further scrutiny. This is, however, beyond the scope 
of this paper. Another point of concern is how the reso¬ 
lution (that is, the number of SPH particles) affects the 
peak temperature. This has been investigated in several 
w orks — see, for instance, th e discussions about this issue 
in lAznar-Siguan et al.l ll2013li . and references therein. It is 
typically found that the peak temperatures may differ by 
a factor of 2 when low- and high-resolution simulations are 
compared, but these temperatures do not show an excessive 
dependence on the adopted resolution when the number of 
SPH particles is larger than ~ 10® — typically the differ¬ 
ences amount to 15%. 


We emphasize that when gdisk = 0.12 is adopted — that 
is, for those mergers for which the eccentricity is high — the 
successive mass transfer episodes only happen when both 
stars are very close. Actually, it turns out that when both 
stars are at closest approach the secondary white dwarf is 
substantially tidally deformed, and the mass transfer episode 
occurs just after the passage through the periastron. How¬ 
ever, soon after this, tidal deformations become smaller, so 
the mass transfer episode stops, and both components of 
the binary system recover their initial spherical symmetry 
for the rest of the orbit. On the contrary, when (/disk = 0.10 
is employed — corresponding to mergers driven by the emis¬ 
sion of gravitational waves — the orbits are almost circular, 
and although the mass transfer episodes take place also at 
periastron, when the secondary is significantly distorted by 


tidal forces, the secondary star remains tidally deformed for 
the entire orbit, and thus the merger process proceeds in a 
smoother way. 

In all four cases the successive mass transfer episodes 
modify the mass ratio of both components of the binary sys¬ 
tem, and this in turn changes the respective orbits. Fig. 0] 
displays the trajectories of the center of mass of both compo¬ 
nents of the binary system. The left panels of this figure de¬ 
pict the trajectories of the centers of mass of the white dwarf 
components. The blue lines correspond to the cases in which 
the core of the AGB star is cold (labelled as CC), while the 
red ones show the cases in which we adopt a high temper¬ 
ature for the core of the AGB star (labelled as HC). In the 
right panels the trajectories of the center of mass of the AGB 
cores are displayed, using the same color coding. The top 
panels of this figure display the case in which gaisk = 0.12. 
As can be seen, in this case the orbits are eccentric, while in 
the bottom panels, corresponding to qaisk = 0.10, the orbits 
are initially nearly circular. As mentioned earlier, for both 
cases it turns out that during the first periastron passage 
the white dwarf transfers mass to the heavier AGB core and 
its radius increases. This radius increase starts a series of 
successive mass transfer episodes that ultimately lead to a 
merger. 

We stress that the number of mass transfer episodes de¬ 
pends not only on the physical characteristics of the binary 
system but also on the number of SPH particles employed in 
the simulations, that is, on the spatial resolution. Although 
our resolution is enough for studying the overall properties 
of the dynamical interaction and the gross properties of the 
merged remnants, it is clear as well that our simulations 
provide only a lower limit to the number of mass transfer 
episodes. For this reason we do not quote the number of 
mass transfer episodes. Note as well that for the case in 
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Figure 6. The evolution of the periods (left panels) and of the eccentricities (right panels) during the merger process. The top panels 
represent the case in which guisk = 0.12 is considered, and the bottom ones those in which guisk = 0.10 is adopted. 


which (/disk = 0.10 it is quite apparent that the last orbits 
are not circular, but instead during the final stages of the 
merger the orbit becomes quite eccentric. This final infall 
phase also takes place in the case of eccentric mergers. How¬ 
ever, due to the difference of scales of the axis, it cannot 
be well seen in the top panels of Fig. [J] Finally, it is also 
worth noting that the trajectories of the centers of mass for 
the case in which cold AGB cores are considered are only 
slightly different of those obtained in the case of hot AGB 
cores. Hence, the effects of the temperature of the AGB core 
on the overall dynamical evolution are very minor. 

Figure [S] shows the mass accreted by the core of the 
AGB star (Mace), in units of the initial white dwarf mass 
(Mwd), as a function of time in units of the initial period of 


the binary system, Tq. As in Fig.jT] the blue lines correspond 
to the case in which a cold merger is considered, while the 
red ones show the evolution for the hot merger. Again, the 
left panel displays the evolution for gaisk = 0.12, whereas 
the right panel shows the same for guisk = 0.10. The mass 
accreted by the core of the AGB star corresponds to the 
total mass of those SPH particles originally belonging to the 
disrupted secondary for which the gravitational attraction 
of the core of the AGB star is larger than that of the white 
dwarf. 

Inspecting Fig. [5] we find that the mass accreted dur¬ 
ing the first mass transfer episode does not depend on the 
temperature of the AGB core, but depends critically on 
the adopted eccentricity. It is also interesting to point out 


© 0000 RAS, MNRAS 000, 000-000 
























10 G. Aznar-Sigudn et al. 




t/To t/T„ 


Figure 7. Average temperature of the core of the AGB star during the merger, for the case in which Qdisk = 0.12 is considered (left 
panel), and for that in which quisk = 0.10 (right panel). 


that for the case in which a hot AGB core is adopted, 
mass is transferred faster than during the first mass trans¬ 
fer episode in successive passages through the periastron, 
for both Qdisk = 0.12 and 0.10. This is more noticeable for 
9 disk = 0.12 — i.e., a merger in which the eccentricity is high. 
Also, it is important to realize that for the cases in which a 
hot AGB core is adopted the amount of mass transferred in 
each one of the episodes is always larger than for the cases 
in which a cold AGB core is considered. All this stems from 
the fact that for a hot AGB core the degeneracy of the ma¬ 
terial of the outer layers is smaller than that of a cold core, 
and thus the core of the AGB star can accommodate more 
accreted mass. 

In all the four cases studied here the mass transferred 
during the second mass transfer episode is larger than that 
transferred during the first periastron passage. This is a con¬ 
sequence of the substantial change of the orbital period af¬ 
ter the first periastron passage, especially for qdisk = 0.12. 
Specifically, for this case the period decreases by « 50% as 
seen in Fig. [O] This is less evident for the case in which 
Qdisk = 0.10 is adopted because in this case the period 
decreases by a modest « 0.2%. All this translates into a 
smoother evolution, also for the successive mass transfer 
episodes, for the gdisk = 0.10 cases. Thus, it turns out that 
the duration of the mergers is shorter for the cases in which 
a hot AGB is adopted, this feature being more pronounced 
for the case in which eccentric orbits are considered. Finally, 
note as well that in each accretion episode there is a small 
decline in the accreted mass. These declines correspond to 
matter, that although being initially accreted by the core of 
the AGB star, is bounced back to the debris region shortly 
after. 

To better illustrate the dynamical evolution of the 
mergers Fig. [6] displays the evolution of the periods (left 
panels) and eccentricities (right panels) of the binary sys¬ 


tems studied here. As before, the top panels depict the evo¬ 
lution of the orbital parameters for gaisk = 0.12, whilst the 
bottom ones show that for the case in which gdisk = 0.10 
is adopted. As can be clearly seen, in all four simulations 
after each mass transfer episode the periods of the binary 
systems decrease and the orbits become circularized. These 
general trends become more accentuated as the evolution 
proceeds. Although for gdisk = 0.10 initially the period de¬ 
creases smoothly in an almost linearly way with a shallow 
slope, after t ~ 35To the white dwarf almost plunges on top 
of the core of the AGB star. This is in contrast with what 
occurs for gdisk = 0.12, for which the decrease in the orbital 
periods occurs faster during the initial stages of the dynam¬ 
ical evolution, in marked steps, and the white dwarf merges 
with the core of the AGB star in very few orbital periods, 
~ 2To. The right panels of this figure show that the eccen¬ 
tricity also decreases for increasing times. For gdisk = 0.12 
this occurs in marked steps, as it occurs with the orbital 
periods, a consequence of the succesive passages through 
the periastron, while for gdisk = 0.10 the eccentricity also 
decreases but oscillates arroung a mean (decreasing) value. 
This is caused by tidal forces. 

Figure [3 shows the average temperature of the core for 
the different cases considered here. To compute this temper¬ 
ature we averaged the temperature of those particles which 
were originally used to model the core, and thus we did not 
take into account the accreted material from the disrupted 
white dwarf. This is an important detail, since the maximum 
temperatures during the evolution of the simulated cases of 
coalescence are not reached in the core of the merged rem¬ 
nant but in the hot corona formed during the interaction 
— see below. As before, the left panel and right panels rep¬ 
resent, respectively, the cases for which gdisk = 0.12 and 
gdisk = 0.10 are adopted. As can be seen, since for eccentric 
orbits the amount of mass transferred during each dynam- 
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Table 1. Some relevant characteristics of the merged remnants. By columns: (1) value of quisk from equation[2 (2) type of run, hot (HC) 
or cold (CC) core of the AGB star, (3) maximum temperature attained during the dynamical interaction, (4) nuclear energy released, 
(5) duration of the merger, (6) Mass of the merged remnant at the end of the simulations, (7) mass of Keplerian disk at the end of 
the simulations, (8) mas of the extended shroud at the end of the simulations, (9) mass of fail-back material, (10) ejected mass, (11) 
maximum temperature of the merged remnant, (12) maximaum angular speed of the hot corona at the end of the simulations. 


Qdisk 

Run 

^peak 

(K) 

-^nuc 

(erg) 

At 

(s) 

Afmr 

(Mq) 

Mdisk 

(Mq) 

-^shroud 

(Mq) 

Mfb 

(Mq) 

Mej 

(Mq) 

Tmax 

(K) 

^max 

(s-l) 

0.12 

HC 

9.23x10® 

1.13x10®® 

18343 

0.88 

0.36 

9.08x10“^ 

2.19x10-2 

1.74x10-2 

3.56x10® 

0.23 


CC 

8.55x10® 

3.54x10®’^ 

24784 

0.88 

0.35 

9.84x10-2 

2 .02x10-2 

1.80x10-2 

3.53x10® 

0.25 

0.10 

HC 

8.74x10® 

1 .12x10®® 

2599 

0.93 

0.39 

4.86x10“2 

2.43x10-® 

7.37xl0-'‘ 

2.97x10® 

0.30 


CC 

8.18x10® 

4.77x10®'^ 

2710 

0.91 

0.41 

4.59x10-2 

2.56x10-® 

8.70x10-^^ 

3.10x10® 

0.32 


ical episode is larger, the temperatures are increased more 
noticeably. Actually, the temperature is not raised homoge¬ 
nously in the entire core. Specifically, the temperature of the 
outer layers of the core is notably increased after every mass 
transfer episode, while the temperature of the very deep in¬ 
terior of the core remains nearly constant. As it should be 
expected, the final temperatures of the hot cores are much 
larger than their cold counterparts. 

Table[T]shows some relevant characteristics of the stellar 
interactions studied in this paper. In columns three to five 
of this table we list, respectively, the maximum temperature 
reached during the dynamical interaction, Tpeak, the total 
nuclear energy released, Enuc, and the total duration of the 
merger. At. The maximum temperature always is reached at 
the base of the accreted layer. The nuclear energy is essen¬ 
tially released by carbon burning reactions and the a chains. 
Note that the nuclear energy released in those mergers in 
which the core of the AGB is hot is typically 25 times larger 
than that released when cold cores are involved. The reason 
for this is that for mergers involving hot AGB cores nuclear 
reactions occur in a region which is significantly larger than 
that in which cold cores are considered. Most of the nuclear 
reactions occur on top of the almost rigid surface of the core 
of the AGB star, as the material of the disrupted white dwarf 
is accreted, compressed and heated. A discussion of the de¬ 
tails of the distribution of the chemical composition in the 
remnant is out of the scope of this paper, and will be pre¬ 
sented elsewhere. However, we mention here that these re¬ 
sults are in line with o ur previous simulations of merging and 
colliding white dwarfs IILoren- Aguilar. Isern fc Garcfa-I^rrol 


coiliQing wUite dwarts llboren- Aguilar, 
I 2 OO 9 I. I 2 OIOI: lAznar-Sieuan etThlboii 


2014), which in turn 


agree with the calculations of Raskin et al.l ||2014| ) and 
iDan et al] ll2014l) . 

Again, we emphasize that the duration of the merger 
not only depends on the physical characteristics of the coa¬ 
lescing stars, but also on the adopted resolution. The reason 
for this is that, as it occurs with DD mergers, the durations 
of the mass transfer episodes depend on the extent to which 
the outer layers of the donor star can be resolved. In most 
practical situations the particle masses are such that the re¬ 
solvable mass transfer rate is al ready a substanti al fraction 
of the total system mass — see iDan et ^ (l201ll ) for a de¬ 
tailed discussion of this issue. Hence, these durations should 
be regarded as indicative, and can be only used to compare 


the different cases studied here, which were computed with 
the same spatial resolution and input physics. 

Figs. H 0 and [S] reveal that the peak temperature is 
reached at t = I.STq and t = 2.4:To for the hot and cold 
cores of the AGB star when gaisk = 0.12 and at t = SS.ITq 
and t = 40.3To for the case in which quisk = 0.10. When this 
occurs the fraction of the secondary white dwarf that has 
been accreted onto the core of the AGB star is > 20% for 
9 disk = 0.12, whereas for qaisk = 0.10 it is < 10%. All this 
has a direct effect on the thermonuclear reaction rates, and 
consequently on the total nuclear energy released during the 
interactions. In particular, while for the case of hot cores the 
total nuclear energy released is « 10^® erg, for the cold ones 
it is much smaller, « 10^^ erg. The nuclear energy released is 
similar for the eccentric orbits and the non-eccentric ones for 
a fixed temperature of the AGB core, despite the difference 
in the peak temperatures achieved during the interaction. 

Finally, we note that eccentric orbits (qaisk = 0.12) re¬ 
sult in mergers with durations of a few hours, much longer 
than those of mergers with qaisk = 0.10, for which the dura¬ 
tions are smaller than ~ 1 hour. However, we remind that 
for the case of eccentric orbits the coalescence occurs in just 
a couple of initial orbital periods, while for circular ones the 
merger lasts for several initial orbital periods. 


4.2 The remnants of the interaction 

In this section we analyze with the help of Figs.l8ll9l ll0l and 
Table[T]the structure of the remnants of the interaction. The 
general appearance of the remnants resulting from the inter¬ 
action is similar for all the cases studied in this paper, and 
moreover is similar to that found in previous simulations 
of the DD scenario for similar masses — s ee, for instance, 
iLoren- Aguilar. Isern fc Garcfa-Berrd (l2009l) . and references 
therein. The remnants of the interactions when simulations 
were stopped consist of a central degenerate object which 
contains all the mass of the core of the AGB star and some 
accreted mass. This central object spins as a rigid body and 
its angular momentum arises from the conversion of the or¬ 
bital angular momentum of the binary system to rotational 
one. This core is surrounded by a corona, which contains a 
sizable fraction, ~ 18% for eccentric mergers and « 27% for 
circular ones, of the mass of the disrupted secondary (the 
white dwarf). The material of the corona is hot — a con- 
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Figure 8. Rotational velocity profiles as a function of the spherical radius in solar units, at times I/Tq = 1.8115 and 2.4476 for the hot 
(solid red lines) and cold (dashed blue lines) AGB cores of the simulation in which guisk = 0.12, and t/To = 40.2378 and 41.9690 for the 
case in which guisk = 0.10 is adopted, respectively. These times correspond to our last computed models. In the left panel we represent 
the profiles for gjisk = 0.12, while the right one corresponds to guisk = 0.10. Note that for large distances the spherical radius is not a 
valid representation of the merged configuration, as the assumption of spherical symmetry is only fulfilled to a good approximation in 
the innermost regions of the merged remnant. The leftmost solid vertical lines show the outer boundaries of the respective degenerate 
remnants. The middle vertical lines show the location of the outer edges of the corresponding coronae. The rightmost vertical lines mark 
the outer edges of the respective Keplerian disks. 


sequence of the material of the disrupted secondary being 
compressed on top of the almost rigid surface of the core of 
the AGB star — and rotates differentially — a consequence 
of the several mass transfer episodes occurring in the dy¬ 
namical interaction. Most of the material of the secondary 
82% and « 73%, respectively, for eccentric and circular 
mergers) that is not accreted on the central object goes to 
form a thick accretion disk, with a Keplerian velocity profile. 
Additionally, there is an extended shroud made of particles 
that have orbits with large inclinations with respect to the 
orbital plane. Also, some SPH particles have highly eccen¬ 
tric orbits and will ultimately fall back onto the primary. 
Finally, we find that very little mass is ejected from the sys¬ 
tem. All these regions are clearly shown in Fig.[8l where we 
plot the angular velocity of the remnant as a function of the 
spherical radius for the remnants of the interactions studied 
here. The vertical lines show from left to right, respectively, 
the locations of the outer edge of the remnant, of the hot, 
rapidly rotating corona and of the Keplerian disk. 

Table [T] lists also some of the characteristics of the 
merged remnant. In column 6 we list the mass of the central 
merged object (Mmr) — that is, the sum of the mass of the 
undisturbed AGB core and of the hot corona — while in col¬ 
umn 7 the mass of the Keplerian disk (Maisk) is listed. The 
precise location of the inner edge of the hot corona is defined 
as the (spherical) mass coordinate for which the angular ve¬ 
locity is no longer constant, while the outer boundary of this 
region is defined as the point where the profile of angular 
velocities of the remnants is the Keplerian. In this table we 
also show the mass of the extended shroud (Mshroud), that 


of the fallback material (Mfb), as well as the mass ejected 
from the system (Mej) — columns 8, 9 and 10, respectively. 
As mentioned, the shroud is composed of material com¬ 
ing from the disrupted secondary that has approximately 
a spherical distribution. The fallback material is composed 
of particles belonging to the debris region that have very 
eccentric orbits as a result of the extreme conditions found 
during the very first stages of the dynamical interaction, 
but are still gravitationally bound to the merged remnant. 
Thus, it is expected that this material will ultim ately fall 
back onto the central remnant at some later time (iRosswod 
120071) . To compute the amount of mass in the debris region 
that will eventually fall back on the c entral object we fol¬ 
lowed the same procedure described in lAznar-Sieuan et al.l 
(|2014|) . Note that there is as well some material that is not 
gravitationally bound to the remnant and might escape, but 
the merger episodes are essentially conservative. We would 
like to remark that we do not look for a possible formation 
of jets and disk winds which, in the DD scenario influece 
the app earance of an explosion if it occurs sho rtly after the 
merger llLevanon. Soker fc Garcfa-Berrd l2015l ). Finally, for 
the sake of completeness we also list the maximum temper¬ 
ature (Tmax) and the maximum angular velocities (wmax) of 
the hot corona, columns 11 and 12, respectively. 

4-2.1 The central compact object 

As can be seen in Table [T] the angular velocities of the 
remnants do not depend appreciably on the adopted tem¬ 
perature of the core of the AGB star, although there are 
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Figure 9. Density (thin lines) and temperature (thick lines) profiles of the central part of the remnants, at the same times shown in 
Fig. El The left panel shows the profiles of the hot (solid red lines) and cold (dashed blue lines) AGB cores of the simulation in which 
9disk = 0.12, while the right panel shows the same quantities for gjisk = 0.10. As before, the abscissa is the spherical radius in solar 
units, and the vertical lines show the same locations explained in Fig. El 


significant differences depending on the value of (/disk. The 
reason for this is that in the coalescence studied here the 
rigid rotation of the core of the AGB star arises from the 
conservation of angular momentum. That is, the orbital an¬ 
gular momentum of the pair is basically invested in spinning 
up the primary and the previously described corona, and to 
form the Keplerian disk. Nevertheless, the angular velocity 
of the remnants of the mergers in which a cold AGB core is 
involved are somewhat larger in both cases. This is due to 
the slightly smaller radii of these configurations (see Fig. [9] 
below), which in turn is a consequence of their larger degen¬ 
eracies. 

We now compare the angular velocities obtained for dif¬ 
ferent values of (/disk. As shown in Table[T] when (/disk = 0.10 
is adopted the angular velocities are larger. In particular, we 
find that for this case the inner regions of the merged rem¬ 
nant rotate uniformly at a speed oj ~ 0.22 s“^, whereas in 
the corona a maximum angular velocity of (Umax — 0.31 s“^ 
is reached. For the case in which (/disk = 0.12 these angu¬ 
lar velocities are, respectively, a; ~ 0.18 s“^, and cumax — 
0.24 s“^, somewhat smaller than those of the case in which 
(/disk = 0.10 is employed. The orbital angular momentum of 
the binary systems studied here is J = 5.96 x 10®° erg s for 
gdisk = 0.12, and J = 4.63 x 10®° erg s when (/disk = 0.10. 
Since the merger is almost conservative (very few particles 
acquire energies large enough to escape from the merged 
remnant) this angular momentum is distributed essentially 
in spinning up the core of the AGB star, in the differential 
rotation of the hot corona and in Keplerian disk. 

For a hot merger their respective angular momenta are 
Jcore ~ 1.4 X 10"^° erg s, Jcorona — 4.4 X lO'^® erg s and Jdisk — 
2.8 X 10®° erg s for gdisk = 0.12, and Jcore — 1.5 x lO'^® erg s, 
Jcorona — 6.7 x lO"'^ erg s and Jdisk — 3.2 x 10®° erg s for 


gdisk = 0.10, respectively. For the case of a cold merger these 
values do not differ much. Hence, the orbital angular mo¬ 
mentum of the eccentric merger is larger than that of the 
circular one, and moreover the angular momenta are dis¬ 
tributed in a different way in both cases. Specifically, the 
angular momentum of the corona for gdisk = 0.12 is notice¬ 
ably smaller than that of the merger in which gdisk = 0.10 
is adopted, but the mass of their respective coronae do not 
differ that much. Also the angular momentum of the newly 
formed disk of the eccentric merger is significantly smaller 
than that of the circular one. Gonsequently, more angular 
momentum is stored in these (outer) regions for the case in 
which (/disk = 0.10 is adopted when compared to the case 
in which Qdisk = 0.12. Finally, for the the simulations in 
which (/disk = 0.12 is employed the ejected mass is consid¬ 
erably larger than that of the case in which (/disk = 0.10 is 
used. This mass carries some angular momentum. All this 
explains why the merged remnant spins at a slower rate for 
gdisk = 0.12. All this is a consequence of the dynamical evo¬ 
lution during the merger. We remind that for Qdisk = 0.12 
the evolution proceeds in abrupt steps, while for qdisk = 0.10 
mass is accreted in a more gentle way. 

Figure |9] shows the density (dashed lines) and temper¬ 
ature (solid lines) profiles, for the case in which qdisk = 0.12 
(left panel) and the for the case in which (/disk = 0.10 (right 
panel). The maximum temperature (Tmax) in all the rem¬ 
nants is very similar, close to 3 x 10® K, being slightly higher 
for the remnants of the eccentric mergers. This can be easily 
understood in terms of the evolution during the last merger 
episode. In particular, it has been already mentioned that 
for eccentric mergers the mass transfer episodes are more 
dramatic than for the case of nearly circular orbits. Hence, 
the matter of the disrupted secondary white dwarf acquires 
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larger accelerations and, consequently, is compressed more 
violently on the surface of the primary star. Thus, the tem¬ 
peratures of the external layers of the merged configuration 
are larger, and a hotter corona is obtained. 

The stronger interaction of the accreted matter with the 
primary during the successive mass transfer episodes leaves 
also clear imprints on the structure of the internal regions 
of the merged remnant. Specifically, for the case of eccentric 
mergers in which a cold core of the AGB star is considered, 
the temperature of the outermost regions of the initially 
isothermal core begins to increase at earlier times during 
the coalescence, when compared to those in which a circular 
orbit is involved. Consequently, by the end of the merger 
process the core of this remnant is no longer isothermal. 
On the contrary, the core in the case of a circular orbit 
remains nearly isothermal — see Fig. |9] Finally, owing to 
the substantial degeneracy of the material of the core of 
the AGB star, the density profiles obtained in all the four 
simulations presented here are very similar. 

4-2.2 The debris region 

It has been already discussed that the temperature of the 
AGB core barely plays a role in the dynamical evolution of 
the merger process, and in determining the structure of the 
merged remnant. Instead, the most important parameter is 
the eccentricity of the initial orbits of both components of 
the binary system. This is true as well for the debris re¬ 
gion formed around the merged remnants. The masses of 
the fallback and of the ejected material are, respectively, 
Mfb — 2.1 X 10“^ Mq and Mej ~ 1.8 x 10“^ Mq for eccen¬ 
tric mergers, while for circular orbits the masses are much 
smaller, Mfb — 2.5 x 10~® Mq and Mej ~ 8.1 x 10“^ Mq. 
The larger fallback and ejected masses in the case of eccen¬ 
tric mergers are a consequence of the longer durations (or, 
equivalently, larger orbital separations) of the mass transfer 
episodes — see Figs. [5]and[6l 

As well, the masses of the respective shrouds are rela¬ 
tively different, Mghroud — 9.2 x 10“^ Mq for guisk = 0.12 
and Mghroud — 4.7 x 10“^ Mq for gaisk = 0.10. This is 
again a consequence of the very different number of mass 
transfer episodes occurring in both simulations. Actually, 
for Qdisk = 0.12 we find that some mass of the disrupted 
white dwarf flows around the region located in the opposite 
direction of the line connecting both stars of the system, 
and acquires high velocities when colliding with the primary 
star, forming this extended shroud. The size of this extended 
region does not depend much on the temperature adopted 
for the core of the AGB star. This can be seen in Fig. 1101 
where the density (left panel) and temperature (right panel) 
contours of the merged remnant in the meridional plane are 
displayed. In this figure the sub-panel with negative values 
of the ^-coordinate corresponds to the case gaisk = 0.12, 
whereas that with positive values shows the same contours 
for the case in which (/disk = 0.10 is adopted. 

As can be seen in Fig. 1101 both the density and the 
temperature decrease as the distance to the central compact 
object increases, independently if we consider the remnants 
of eccentric or circular mergers. However the decline rate of 
both the density and the temperature for eccentric mergers 
is smaller than for circular orbits, resulting in a more ex¬ 
tended debris region, which moreover is hotter. Additionally, 


the density and temperature profiles also have a different 
behavior near the polar regions of the central remnant. For 
mergers arising from circular orbits there is a region near the 
z-axis (the vertical axis in this figure) where the density and 
the temperature contour lines converge towards the poles, in 
contrast with what happens for the case of eccentric merg¬ 
ers, for which this region is absent, and the contour lines are 
perpendicular to this axis. These differences between both 
profiles are a consequence of how the material of the dis¬ 
rupted secondary white dwarf is distributed on top of the 
primary (the core of the AGB star) during the dynamical 
interaction. 

The final infall of the disrupted secondary onto the un¬ 
altered primary for guisk = 0.10 seen in Figs. [3] and [4] does 
not allow matter to be redistributed on the primary, but in¬ 
stead concentrates it in the orbital plane. Besides, the circu¬ 
lar density contours that appear at the right of the remnant 
correspond to a spiral arm of matter infalling to the primary. 
Since SPH codes are not able to follow properly angular 
momentum transport, we stopped the simulations at a rea¬ 
sonable time after the complete disruption of the secondary. 
In particular, for gaisk = 0.12 we stopped the simulations 
at I/Tq = 1.8115 and t/To = 2.4476 for the hot and cold 
merger, respectively. For gaisk = 0.10 these times are, respec¬ 
tively, t/To = 40.2378 and t/To = 41.9690. In this case, the 
secondary had no time to be totally distributed over the pri¬ 
mary and the spiral arm formed during the infalling phase is 
still present. The more flat structure for qaisk = 0.10, which 
is more similar to that usually found in the simulations of 
the DD scenario, might lead to jets and disk winds that 
might alter the observations if explosion occurs s hortly af¬ 
ter merger (ILevanon. Soker fc Garcfa-Berra l2015ll . Finally, 
the outer edges of the Keplerian disks formed during the 
merger process are located where the density distribution 
in the equatorial plane has no longer axial symmetry. In all 
cases this occurs at ~ 0.06 7?q. 


5 SUMMARY AND CONCLUSIONS 

In this paper we have simulated the coalescence of a white 
dwarf with the core of a massive AGB star, in the context 
of the so-called core-degenerate scenario for SNe la. Specifi¬ 
cally, we have computed the merger of an otherwise typical 
white dwarf of mass 0.6 Mq with the core of an AGB star 
of m ass 0.77 Mq, which is a typical case (iKashi fc Sokeil 
This has been done for two temperatures of the AGB 
core, T — 10® K and 10® K, and for two assumptions about 
the mass of the disk formed during the common envelope 
phase. In particular, we have studied a first case in which 
the mass of the disk is Mdisk = 0.12(Mcore + Mwd), which 
corresponds to the critical mass ratio for which a merger 
is guaranteed as a consequence of the interaction with the 
disk formed during the common envelope phase. The second 
case studied here corresponds to a binary system surviving 
the common envelope phase with a mass ratio 0.10, which is 
below the critical one, and for which the merger is driven by 
the emission of gravitational waves. The first of these cases 
results in a pair with a highly eccentric orbit, whereas for 
the second one the orbits of the members of the binary sys¬ 
tem are circular, and similar to those expected for the DD 
scenario for SNe la. 
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Figure 10. Left panel: density contour lines of the remnants across the meridional plane. Right panel: temperature contour lines across 
the same section. These density and temperature maps correspond to our last computed models, and the evolutionary times are the same 
of Fig. 1^ The remnants represented with negative x coordinates correspond to the case in which guisk = 0.12, while for positive values 
of X the case in which gjisk = 0.10 is represented. For both cases the outermost contour of the density profile corresponds to a density 
of 1 g cm“® and the successive contours increase logarithmically inwards in constant steps of logp = 0.62. The contour of maximum 
density corresponds to ~ 3.8 X 10® g cm“®. For the temperature the outermost contour corresponds to ~ 3 X 10® K and the successive 
contours increase inwards logarithmically in steps of logT = 0.33. The maximum temperature contour corresponds to ~ 1.4 X 10® K. 
This temperature is reached in the hot corona, while the innermost regions of the merged remnant in the case of a core with small 
temperature merger are colder — see Fig. [7] See the online edition of the journal for a color version of this figure. 


Our SPH simulations show that temperature of the core 
of the AGB star has little influence on the overall dynamical 
evolution of the merging process and on the characteristics 
of the merger remnant, but the initial eccentricity of the 
binary system has a large impact on both the dynamical 
evolution of the merger and on the properties of the merged 
object. In particular, in both cases the white dwarf is to¬ 
tally disrupted in a series of mass transfer episodes, which 
occur at closest approach. However, the duration of the co¬ 
alescence is notably different in both cases. Specihcally, we 
found that for the case of an eccentric merger the disruption 
of the secondary white dwarf occurs in very few orbital pe¬ 
riods, whereas for the case of the mergers in which circular 
orbits are adopted this process lasts for many orbital peri¬ 
ods. Nevertheless, the initial orbital periods in both cases are 
very different — 10,126 s and 65 s, respectively — so the to¬ 
tal time required to disrupt the secondary white dwarf turns 
out to be larger for the case of eccentric orbits. 

In all the four cases studied in this paper the coa¬ 
lescence process is almost conservative, and little mass is 
ejected from the system, ~ 2 x 10”^ Mq at most. Neverthe¬ 
less, eccentric mergers eject ~ 25 times more mass than 
that ejected by binary systems with circular orbits. The 
peak temperatures achieved during the merging process are 
rather high, of the order of 8 x 10® K, but not enough to 
power a prompt detonation upon merger. As for the final 
temperature of the merger remnant it is interesting to re¬ 
alize that for the cases studied here the temperature pro¬ 
file peaks off-center. The possibility that carbon is ignited 


in the remnant, thus producing a rapidly s pinning oxygen- 
neon w hite dwarf, has been discu ssed bv ISaio Nomotd 
ll2004h , lYoon fc Langerl (|2004| ~I , and lYoon fc Langerl ( 2005ll . 


Our results do not support a prompt off-center ignition, 
since the peak temperature obtained at the end of our sim¬ 
ulations is too low, and moreover the temperature peaks 
in a region where the densities are relatively small. How¬ 
ever, even though carbon in our models is not ignited im¬ 
mediately during the merger process, it may occur at later 
times, when the high entropy material begins to cool and 
compress the material near the temperature peak. Never¬ 
theless, this depends on the not yet well known viscous 
timescale of the disk, which characterizes the transport of 
disk mass inwards and angular momentum outwards, and on 
the thermal timescale of the merger product. Thus, the post- 


(Ivan Kerkwiik. Ghang & Justhamll2010l: Schwab et al.l 

2012 ; 

Shen et al. 

120121 : 

Garcfa-Berro et al. 20121: Ji et al. 

2013; 

Zhu et al.l 1 

20131: iBeloborodovI |2014|l but a clear consensus 


on whether a powerful carbon ignition can develop at later 
times has not emerged yet. 


The structure of the merged remnant consists of a cen¬ 
tral compact object containing all the mass of the core of 
the AGB star and some mass from the disrupted secondary 
star, which is surrounded by a hot, differentially rotating 
corona. This corona is made of about 18% of the mass of 
the disrupted white dwarf, for the case of eccentric merg¬ 
ers, and somewhat larger (~ 27%) for circular ones. Finally, 
surrounding this central compact object there is a heavy. 
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rapidly-rotating Keplerian disk, and an extended shroud of 
material which are formed by the rest of the mass of the 
disrupted white dwarf that has approximately ellipsoidal 
shape. This configuration is essentially the same found for 
DD mergers. However, the structure of the debris region is 
different depending on the initial eccentricity of the merger. 
For mergers with initially circular orbits the matter of the 
disrupted white dwarf concentrates in regions near the or¬ 
bital plane, and little mass is found in the polar regions, 
whereas for eccentric mergers part of the mass of the dis¬ 
rupted white dwarf is distributed along the central compact 
object, and the merged configuration adopts an ellipsoidal 
shape. In the debris region there is also some fraction of SPH 
particles with highly eccentric orbits. This material will ul¬ 
timately interact with matter in the Keplerian disk and will 
fall back onto the central compact remnant. The mass of this 
fallback material is typically or the order of 2 x 10“^ Mq 
for the case of the binary systems with eccentric orbit, and 
ten times smaller for those with circular ones. In all cases, 
this central object rotates as a rigid body, with typical an¬ 
gular velocities ~ 0.18 s“^ for the case of eccentric mergers 
and ~ 0.22 s~^ for circular ones. These rotation velocities 
arise from the conversion of orbital angular momentum to 
rotational one, as it happens for DD mergers. 


We would like to emphasize that if the disk, fallback 
material, and the material of the extended shroud are 
hnally accreted, as it is commonly thought, the mass of 
the remnant would be « 1.35 Mq. If this is the case, after 
cooling and losing angular momentum this remnant might 
eventually explode. Of course, assessing such possibility 
requires modelling the viscous phase of the merger, which 
transports disk mass inwards and angular momentum 
outwards. This phase has typical timescales ranging from 
10 ^ s to 10 ^ s, during which the disk might lau nch jets 
and winds (|Levanon. Soker fc Garcfa-Berrol l2015l l , and it 
cannot be followed using SPH techniques. If, nonetheless, 
all this mass is accreted and the remnant explodes the delay 
would be very long, up to millions of years, or even longer 
Jllkov fc Sokerj[2012l l. In particular, for accretion of helium 
onto a carbon-o xygen white dwarf in a stable reg i me, in 
a recent study IPiersanti. Tornambe fc YungelsonI (|2014l 'l 
found a delay of up to millions of ye ars, in agreement with 
previous calcula ti ons o f t his kind dNom oto &: IbenI 19851: 
Saio fc Nomotol l2004l : lYoon. Podsiadlowski fc Rosswod 


200?! : IShen fc BildsterJ Ibood l . It remains to be checked if 

the same might occur for the CD merger process. 


In summary, this first set of simulations of the CD sce¬ 
nario for SNe la demonstrates the need of exploring in more 
detail the effects of the eccentricity (and possibly of the tem¬ 
perature) on the overall dynamical evolution of the merger 
and on the properties of the merged configurations obtained 
in this scenario, especially when more massive stars are in¬ 
volved. The simulations presented here are a first important 
step in this direction. Future works will undoubtely allow 
us to characterize and place tight constraints on the most 
salient features of this scenario, and may possibly help in 
quantifying the fraction of these mergers that contributes to 
the total SNe la rate, thus filling a gap and complementing 
what we already learned from the existing simulations of the 
DD scenario. 
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